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Electron trapping by electric field reversal in glow discharges 
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The phenomena of electric field reversal in glow discharges is discussed. Several models are 
described and the link between field reversal and the discharge structure is analyzed. 
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length. They obtained as well the fraction of ions ar- 
riving at the cathode and the magnitude of the plasma 
maximum density. 

I INTRODUCTION Technological application of gas discharges, particu- 

larly to plasma display panels, needs a better knowledge 
of the processes involved. 

The phenomena of field reversal of the axial electric 
field in the negative glow of a dc discharge has been dis- 
cussed since a long time in literature P, [^. In a founding 
paper Druvestyn and Penning [H hypothesized its exis- 
tence. Lawler et al. 3 assumed its existence when evalu- 
ating the current balance at the surface of a cold cathode 
discharge and estimated the position of the field reversal 
by linearly extrapolating the optical data of atomic tran- 
sitions. The determination of the position where the field 
reversal occurs is of great importance since the fraction 
of ions returning to the cathode depends on its existence 
and location, and it is necessary for correctly determine 
the conditions for a self-sustained discharge. 

Kolobov and Tsendin Q have shown that the first 
field reversal is located near the end of the negative glow 
(NG), near the position (although located slightly to the 
cathode side) where ions density attains its greatest mag- 
nitude 6]. If the discharge length has enough extension 
and the pressure decrease to lower values it appears a sec- 
ond field reversal on the boundary between the Faraday 
dark space and the positive column (PC) 0, Q. 

Moreover, Kolobov and Tsendin explained how ions 
produced to the left of the first reversal location move 
to the cathode by ambipolar diffusion - helping to main- 
tain the glow by secondary electron emission - and ions 
generated to the right of this location drift to the anode. 

More recent work [tI, [§[ presenting a comparison of 
experimental data and the predictions of a hybrid fluid- 
Monte Carlo model also supports the view that the point 
where the field is extrapolated to zero is practically co- 
incident with the maximum of the emission (even when 
j/p^ scaling is no longer valid). Those characteristics 
were experimentally observed by laser optogalvanic spec- 
troscopy Q . For a detailed review see also . 

Boeuf and Pitchford [l^ with a simple fluid model 
gave an analytical expression of the field reversal loca- 
tion showing its dependency solely on the cathode sheath 
length, the gap length, and the ionization relaxation 



'Electronic address: ' mpinheiro®ist.utl.pt| 



II. STRUCTURE OF A GLOW DISCHARGE 

The phenomenology of DC glow discharges are well 
described, for example, by Roth [Till, where the existence 
of the field reversal in the negative glow (NG) is already 
clearly shown. To be complete, we will recall now the 
most significant features with interest to our problem. 
In the regime of a normal glow discharge the voltage and 
the current density are both practically independent of 
the total current fiowing in the discharge tube. 

In a crude picture, energy is continuously transferred 
from a high voltage DC power source to the electrons 
originating from the cathode that accelerate absorbing 
energy from the field eE, ionizing, exciting and undergo- 
ing elastic collisions with heavy particles and other elec- 
trons. Electrons disappear basically by recombination 
and diffusion to the walls. Secondary electrons from the 
cathode provides an additional source of electrons and 
they have a prominent role in the discharge maintenance. 
The electric field decreases in the sheath, vary slowly in 
the plasma and increases again in the anode sheath, al- 
though not attaining such high value as in the cathode. 
All plasmas are separated from the walls (both conduct- 
ing or non-conducting) by a sheath. The higher ener- 
gies are attained by electrons which bombard far more 
frequently the walls than ions. The slower ions remain- 
ing behind build-up a positive plasma potential relative 
to the wall with a magnitude of the order of the elec- 
tron kinetic temperature. In consequence ions are eas- 
ily drawn off through the sheath by an external circuit 
as a current / while electrons are repelled by the walls 
(and electrodes), only escaping the most energetic ones. 
As we proceed from the cathode to the anode we recog- 
nize several regions, first observed by Michael Faraday in 
the 1830's, with interest to the present study (see also 
Refs. ^m)- 

• Aston dark space - The electric field is very high 
and primary and secondary electrons outnumber 
ions. However, electrons have not still attained 
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enough energy to excite neutrals and they have low 
density making this region to appear dark; 

• cathode glow - electrons have already enough en- 
ergy to excite neutrals and ions density increases; 

• cathode (Crookes, Hittorf) dark space - Electrons 
are either too fast or too slow to excite the gas. 
The probability for fast electrons to excite are not 
negligible, however, and that's why this regions is 
slightly brighter than the Aston dark space. In 
this region electrons have enough energy to ionize 
the gas and the ions production rate attains high 
magnitude, the dominant charge being positive. It 
is in this region that the primary electrons loose 
a great fraction of its energy in ionizing collisions 
and the generated secondary electrons appear with 
lower energy. Slow electrons are produced in this 
region. A moderate electric field dominates; 

• cathode region (CR) - This is a transition region 
between the cathode dark space and the negative 
glow where most of the voltage drop (cathode fall) 
of (pc 'Ai X = dc- The axial length of the cathode 
region dc adjust its value in order to maintain a 
minimum value of the product dcP, the Paschen 
minimum; 

• negative glow (NG) - The electric field is here very 
low, but its luminosity results from the electrons 
that have been accelerated in the CR and ionize 
and excite neutrals in the NG (that's why the lu- 
minosity is more intense on the CR side). As elec- 
trons lose their energy, recombinations processes 
are more probable, mainly in poliatomic electroneg- 
ative gases to which electrons attach more easily. 
Electrons carry almost all the entire electric cur- 
rent in this region. Ionization in the glow can con- 
tribute significantly to the total ion current in the 
glow due to backscattered ions from the NG enter- 
ing the cathode fall [l^ . 

• Faraday dark space (FDS) - This region is located 
immediately to the right of the NG. Electrons en- 
ergy is low and their density decreases due to re- 
combination and radial diffusion. However, in gas 
discharge tubes with a radius much smaller than 
the discharge length, the electric field start to grow 
again and electrons gain energy giving birth to the 
positive column. 

III. EXPERIMENTAL OBSERVATION OF 
FIELD REVERSAL 

The phenomenon of field reversal was experimen- 
tally observed by laser optogalvanic spectroscopy (LOG) 
firstly by Gottscho et al. [6]. The_LOG technique used 
was reported by Walkup et al. [l4| and is based on a 
change of ion mobility on excitation by a dye laser. This 



means that an excited ion state with a larger (smaller) 
mobility will induce a transient current increase (de- 
crease) on excitation. The sign of the LOG signal de- 
pends only on the sign of electric field reversal. Gottscho 
et al. (6j] at the same time recorded the N2 {X^'E,'^) ions 
excited from its lowest and first excited vibrational levels 
(exciting the B'^E+ - X'^T,,g+ (0,0) and (0,1) bands, re- 
spectively) using laser-induced fluorescence (LIF) along 
the plane parallel electrodes of their reactor. They noted 
that at low pressure the LOG signal changes sign near 
the position of the iV^ maximum while being located 
slightly to the cathode side. At higher pressure the field 
reversal occurs toward the anode side of the density 
maximum. 

The LOG and LIF spatial profiles resulting from the 
excitation of the B^J:+ - X'^T..,g+ (0,0) and (0,1) bands 
decrease first strongly to negative values due to the falloff 
in ion density, as it is registered by LIF. This is the pre- 
shcath pj| . a quasi-neutral region between the plasma 
and the sheath characterized by small electric fields gra- 
dients (usually smaller than 1 V/cm). At a given posi- 
tion ions practically remain undetected by LIF, the LOG 
signal diminishes rapidly and this is a signature of the 
beginning of the cathode fall region. 



IV. PHYSICAL MECHANISM OF THE 
ELECTRIC FIELD REVERSALS IN GLOW 
DISCHARGES 

Electric field reversals operate as a self-regulator mech- 
anism of glow discharges. If electrons were in hydrody- 
namic equilibrium with the local electric field and were 
created in the cathode sheath followed by a fast diffu- 
sion to the walls and acceleration to the anode, then ions 
loss would occur through migration to the cathode and it 
will be no need for field reversal. This situation would be 
similar if the ionization main channels were ion and fast 
neutral impact ionization instead of electron impact ion- 
ization. However, electrons are not in equilibrium with 
the local electric field. This situation results from the 
complex structure embodied in the sheath, pre-sheath 
and plasma regions (see Ref. [IH for details). The electric 
field decreases largely from the cathode surface within an 
adjacent region where the charged-particle density is low 
(e.g., [l5|). In this region electrons gain energy being 
strongly accelerated toward the bulk plasma attaining 
typically kinetic energies in the range 10 ^ 10"^ eV. Thus, 
they are responsible of the excitation and ionization pro- 
cesses in the quasi-neutral region we call plasma. Clearly 
their energy is determined by the electron-energy spec- 
trum and not by the local value of the electric field . 

Under typical gas discharge conditions it is possible to 
consider two groups of electrons: the group of slow and 
fast electrons. The majority of them belongs to the slow 
group of electrons; they are created in the quasi-neutral 
(bulk) plasma region through ionization and their kinetic 
energy does not exceed significantly the energy of the first 
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excitation level ei The low electric field prevalent 

in the quasi-neutral region with a vast majority of slow 
electrons (with energies in the range lO^'^ ^ 10 eV) is 
favorable to increase electrons loss while at the same time 
reducing ions loss. But as the total current density 



tion for fast electrons F(v,x) is written under the form 



J = Je 



-eneVe 



(1) 



must be constant along the discharge a mechanism of 
field reversal appears in order to enhance ions loss and 
reduce slow electrons loss. At the position where field 
reversal takes place the plasma potential presents a min- 
imum value [6] . This situation creates a kind of two con- 
tiguous reservoirs: to the left of the field reversal, ions 
generated by fast electrons flow to the cathode, where 
hitting the electrode surface induce secondary electrons 
emission at the cathode (and thus helping maintaining 
the glow); to the right of the field reversal, ions are ac- 
celerated toward the anode. As electrons diffuse more 
rapidly to the walls than ions, the walls along the posi- 
tive column and near the anode sheath charge negatively, 
while along the cathode fall they charge positively (see, 
for example, Ref. This charge build-up deform the 

equipotentials surfaces that, instead to remain parallel to 
the electrodes, show the concave side facing the anode. 



KOLOBOV AND TSENDIN NONLOCAL 
KINETIC MODEL 



Due to the complex self-consistent problem that is at 
stake, Kolobov and Tsendin endeavored to develop an 
analytical approach with a formidable com plex ity. Since 
the papers of Emeleus and co-workers [l^ll9j| electrons 
have been divided in primary, secondary (or interme- 
diate) and ultimate (or trapped) and so did Kolobov 
and Tsendin. The highly energetic electron created in 
the cathode sheath of a DC discharge have kinetic en- 
ergies well above the first excitation potential e*. Ion- 
ization and excitation processes are created by them and 
they carry the electron current in the sheath and nearby 
plasma region. The plasma region is mainly populated 
by trapped electrons that do not contribute to the cur- 
rent. In the NG and FDS the current is mainly carried 
by (untrapped) intermediate electrons with energies be- 
low e* . The discharge gap is divided in two regions: the 
sheath and the (quasi-neutral) plasma region. 

The fast electrons are described by a continuous- 
energy-loss model neglecting scattering 



dw 



-NC{w) - eE{x), 



(2) 



where w is the electron kinetic energy, N is the neutral 
particle density, C{w) is the energy-loss function, ^ is the 
fast-electron path along its trajectory and x is the spatial 
coordinate. 

Neglecting scattering, then ^ = x and the kinetic equa- 



dF eE{x) OF d NC{w) 



F = J{x,v). (3) 



dx m dv dv r 

Here, J{x, v) is a source term and it was assumed C{w) ~ 
Cq = const. Its integration is performed at constant 
energy 



w + NCqx — e(l){x) = w — e(j){x), 



(4) 



where 4>{x) is the electrostatic potential obtained for 
linear-electric field profile in the sheath. The fast elec- 
trons start riding on an effective potential, ecfi^x), that 
attains a minimum in a point where the total force is 
zero, eE{d) = NCq (see Fig. 2). For example, for He, 
Co = 1.5 X 10^^^ eV cm^. Solving Eq. [3] it is obtained 
the fast-electron current 

/•OO 

jf{x) = e / vdvF{v,x) = 

jJx) = l ereexp(aa;), ^ < a; < 

\ jf{xo{e ^ e(j){x))), d < X < xi{e4>c), 

where a = NCq/eq, with Eq denoting the constant energy 
loss per ion-electron pair and Fe is the electron flux from 
the cathode. Notice that for x < d the total electron 
current je is equal to the fast-electron current jf and 
that, in particular, we have 



dje 

dx 



(6) 



Here, I(x) is the ionization density. Eqs.[5]and [5] con- 
stitute a generalization of the Townsend approach with 
a constant ionization coefficient a — const. 

At X > xi{e4>c) there is no ionization source since the 
fast electrons do not reach this region; this region limits 
the boundary between the FDS and NG (see Sec. 2 for 
its characterization) and the length of the NG is given 
qualitatively by J? = e^c/o^^o- 

For slow electrons with w < e* the continuous-energy- 
loss is no longer valid and it is essential to take into ac- 
count the essential discrete character of the energy loss, 
where elastic and electron-electron collisions are the only 
energy-relaxation mechanisms. If L is bigger than the 
electron mean free path A, then the electron energy distri- 
bution function (EDF) is isotropic and the two-term ex- 
pansion in spherical harmonics holds. Hence, taking into 
account the following contributions: i) degradation of 
fast electrons injected into the plasma from the cathode 
sheath; ii) slow electrons generation by fast electrons; iii) 
slow electrons superelastic collisions with molecular low 
vibrational levels or others, with kinetic energy slightly 
exceeding e*, the kinetic equation for the isotropic part 
of the EDF gives for the trapped electrons a Maxwell- 
Boltzmann distribution of the type: 
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nm[exp(-e/re) - exp(-e0a/T'e)], 



(7) 
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FIG. 1: a) Total force acting on an electron under constant 
energy-loss; b) effective potential energy of electrons, e4>{x)) 



where <f>a denotes the anode potential and jim is the max- 
imum trapped electron density at the point of field re- 
versal E{x) = 0. The electron temperature is obtained 
by solving the integral energy balance for trapped elec- 
trons, but in a first approach we have e4>f, ^ Tg. For 
example, by laser diagnostics experiments [l^ the slow 
electron temperature and density are = 0.2 eV and 
rim = 1-2 X 10^^ cm'^in He at p = 3.5 Torr, with an 
electrode gap L — 0.62 cm and j — 0.846 mA/ cm^. 
The trapped electrons give the main contribution to the 
plasma density 



ns(x) 



An 



v{e,x)fQ{e,x)de, 



(8) 



while giving zero contribution to the slow-electrons cur- 
rent density: 
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e(p{x 



v{£,x)D{e,x)^^de. (9) 



The js is the the total electron current in the FDS rep- 
resented by the transport of the intermediate untrapped 
electrons. In the low-pressure limit the trapped electrons 
are expected to be isothermal and to be governed by the 
ambipolar diffusion equation: 



d f ^ dn\ ^, , n 
dx \ dx J T 



(10) 



Here, r = A? /Da- Eq.[TO]can be solved using the Green's 
function approach with appropriate boundary conditions. 
It gives the minimum Xm'- 



I xl{x)dx ^ (L + d) [ 



I{x)dx. 



(11) 



The position of Xm determines the ratio of ions return- 
ing to the cathode. If the FDS length (L — R) exceeds 
[R — d) then Xm R and the majority of ions gen- 
erated in the plasma return to the cathode; if, on the 
contrary, R — d > L — R this means that fast electrons 
penetrate deep in the plasma, the majority of them at- 
taining the anode. In their calculation the position of 
maximum ions density coincides with the field reversal 
location. Unfortunately, Kolobov and Tsendin do not 
give explicit expression for Xm, introducing an error in 
Eq. [11] (see remark in Ref. 

To resume the previous work by Kolobov and 
Tsendin ^ , the study of nonlocal phenomena in electron 
kinetics of coUisional gas discharge plasma have shown 
that in the presence of field reversals the bulk electrons in 
the cathode plasma are clearly separated in two groups of 
slow electrons: trapped and free electrons. Trapped elec- 
trons give no contribution to the current but represent 
the majority of the electron population. The first field 
reversal it was shown qualitatively to be located near the 
end of the negative glow (NG) where the plasma density 
attains the greatest magnitude. If the discharge length 
is long enough, it appears a second field reversal on the 
boundary between the Faraday dark space and the pos- 
itive column. Also, it was shown that ions produced to 
the left of this first reversal location move to the cathode 
by ambipolar diffusion and ions generated to the right of 
this location drift to the anode. For a review see also M. 



VI. BOEUF AND PITCHFORD FLUID MODEL 

Boeuf et al [l^ with a simple fluid model gave an an- 
alytical expression of the field reversal location which 
showed to depend solely on the cathode sheath length, 
the gap length, and the ionization relaxation length. 
They obtained as well a simple analytical expression giv- 
ing the fraction of ions returning to the cathode and the 
magnitude of the plasma maximum density. 

Their model assumes a spatial profile of the ioniza- 
tion rate St{x) = Ne{x)i>i{x) in a plasma being formed 
in a parallel plate geometry. Monte Carlo calculations 
show that the ionization rate is linear in the sheath and 
decreases exponentially in the NG [2l| and so they as- 
sumed 



Srix) 





s exp(- 



,for X < dc 
^) ,fora;>4 



(12) 



Here, Ae is the energy relaxation length of the high elec- 
trons energy entering the NG accelerated in the sheath 
while s is a function of the voltage (j^c gas mixture and 
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is proportional to the flux of secondary electrons leaving 
the cathode T^o (here, e (electrons) and i (ions) denotes 
the specie). They assume that slow electrons can be de- 
scribed by a continuity equation and a momentum trans- 
port equation in the drift-diffusion approximation (and 
thus considering drift energy negligible with respect to 
thermal energy): 



dx 



Srix), 



with 



drie 



dx ' 



(13) 



(14) 



where pLe^i and De,i denote the electron and ion mobility 
and diffusion coefScients respectively. Quasi-neutrality 
{n = Tie ^ Ui) is assumed in the NG region. The total 
current density is given by the expression 



Jt = e{-(t)e + 4>i). 



(15) 



From this equation, together with Eq. [HI it is obtained 
the electric field 



E{x) 



Jt 



D, 1 dn ^^g^ 



eUeilie + Ml) tJ'e+ fJ-i n dx ' 

where the electron and ion fluxes are now given by 



and 



(f>eix) 



(f)i{x) 



e(/Xe -I- yUj) 



D — 

dx 



Da 



dn 



(17) 



(18) 



e{iie+fJ'i) dx' 
Here, Da is the ambipolar diffusion coefflcient given by 



Me + Ml 



(19) 



Integrating Eq. ll3l and taking into account the ionization 
source terms in Eq. [12] it is obtained 



re(a;) 



Feo ,for x < dc 

Feo + sAe[l — a(a::)] ,for x > dc 



(20) 



using the simplifying notation a{x) — exp[— (x — dc)/ Xe]- 
At some position x^ it can possibly occur a fleld rever- 
sal, located at a boundary that split the ion flux at the 
right side toward the anode and keeping electrons con- 
flned in the plasma region while at the left side directing 
the ion flux toward the cathode. Hence, a condition can 
be stated of a total current at some yet unspecifled po- 
sition xq (in this Sec. we use this notation for the fleld 
reversal location) being equal to the electron current: 



Jt 



-ere(a;o) = -e{Feo + sXe[l - a{xo)]}- (21) 




FIG. 2: Position of the field reversal, ^ 
the adimensional parameter A. 



as a function of 



From Eq. [ITj and Eq. [20l using Eq. [STJ it is obtained the 
following flrst order differential equation for x > d^. 



l+ — [l-a{x)] 

(PeO 



Me 



s\ , . M Dn dn 

(Me + Mi) V 'PeO 0eO dx 

(22) 

Integrating Eq. [22| between dc and L with appropriate 
boundary conditions, n{dc) = n(d) = 0, under the ap- 
proximation fii <^ ^e, it is obtained 



a(a;o) = 



A 



d — d, 



■[l-a{d)]. 



(23) 



Introducing the dimcnsionless parameter A ~ X/(d~' dc), 
representing the ratio of the relaxation length to the dis- 
tance between the anode and the CR boundary, Eq. [221 
can be written under the form 



Xq 



L - dc 



-Aln[A(l -exp(-l/A)]. 



(24) 



The above equation shows that, by one side, the position 
of the field reversal xq depends on the cathode sheath 
length dc, the gap length L and on the electron energy 
relaxation length A^. By other side, dc and A^ depend on 
the cathode sheath voltage and gas mixture. As can be 
seen in Fig.[2l when A^ <C i (low values of A) the position 
of the field reversal is near the plasma-sheath boundary; 
when Ae ^ L (large values of A) the position of the 
field reversal is displaced to the mid-distance between 
the gap length and the plasma-sheath boundary. This 
last case occurs for an obstructed discharge, i.e, when 
the gap length is less than the Paschen minimum dc at 
the Paschen minimum. 

The fractions of ions to the cathode can be determined 
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0,8 — 




FIG. 3: Fraction of ions returning to the cathode as a function 
of the adimensional parameter A. 



using the relationship 



V = 



lo " ST{x)dx (j)e{xQ) ~ 0eO 
ST{x)dx ~ ML) - (f>eO ■ 



(25) 



We can give an explicit form to this ratio using Eq. [20] 
and Eq. El 



1 



1 - exp(-A) 



A. 



(26) 



This fraction depends only on A and, as Fig. [3] illustrate, 
its values are between 0.5 and 1. 

Fig. [3] shows that when the electron energy relaxation 
length A J is smaller than the gap length, rj ~^ 1, the ma- 
jor fraction of ions return to the cathode. This is true if 
radial diffusion to the walls or recombination is not ac- 
counted for. In the condition of an obstructed discharge 
Ae 3> L, then rj — > 0.5. Taking into account radial dif- 
fusion and recombination will lower the fraction of ions 
returning to the cathode. 

Integrating Eq. EH between x ^ dc and a given x, 
considering appropriate boundary conditions, n{dc) = 
n(L) = 0, assuming ^ /ig and using Eq. E21it is ob- 
tained readily 



n{x) = ^(1 



exp(-l/A)) 



1 - cxp(-x/A) 
1 - exp(-l/A) 



with X 



-Ac. 



. (27) 

jjzTd'- The plasma density attains a maximum 
at the position where occurs the field reversal a;o. This is 
inconsistent with Eq. 1161 however, since here it was ne- 
glected the ions mobility in regard of the electrons mobil- 
ity. In fact, the maximum of the plasma density is located 
near the position of the field reversal, but its maximum 



should be located to the cathode side of the field rever- 
sal. Marie et al. 0, [1] showed a perfect agreement of the 
analytical formula given in Eq. EH with an hybrid Monte 
Carlo-fluid model. Kudryavtsev and Toinova [20*1 deter- 
mined the position of the field reversal from the position 
of maximum ion density in the plasma. Their approach is 
quite simple giving simple analytical formulas in a short 
(without positive column) DC glow discharge making use 
of the ionization rate coefficients for fast electrons in the 
field. 



VII. A DIELECTRIC-LIKE MODEL OF FIELD 
REVERSAL 

In this section we introduce a quite simple dielectric- 
like model of a plasma-sheath system 1231. This approach 
have been addressed by other authors j23l. [2^ to explain 
how the electrical field inversion occurs at the interface 
between the plasma sheath and the beginning of the neg- 
ative glow. The aim of this Letter is to obtain more infor- 
mation about the fundamental properties related to field 
inversion phenomena in the frame of a dielectric model. 
It is obtained a simple analytical dependence of the axial 
location where field reversal occurs in terms of macro- 
scopic parameters. In addition, it is obtained the magni- 
tude of the minimum electric field inside the through, the 
trapped well length, and the trapping time of the slow 
electrons into the well. We emphasize in particular the 
description of the dielectric behavior and do not contem- 
plate plasma chemistry and plasma-surface interactions. 

The analytical results hereby obtained could be useful 
for hybrid fluid-particle models (e.g., Fiala et al. [25j). 
since simple criteria can be applied to accurately remove 
electrons from the simulations. 

On the ground of the stress-energy tensor considera- 
tions it is shown the inherent instability of the field in- 
version sheath. The slow electrons distribution function 
is obtained assuming the Fermi 26] mechanism respon- 
sible for their acceleration from the trapping well. 

Lets consider a plasma formed between two parallel- 
plate electrodes due to an applied dc electric field. We 
assume a planar geometry, but extension to cylindrical 
geometry is straightforward. The applied voltage is Va 
and we assume the cathode fall length is / and the nega- 
tive glow -|- eventually the positive column extends over 
the length Zg, such that the total length \s L = 1 + We 
have 



Va = lEs + loEp 



(28) 



where Eg and Ep are, resp., the electric fields in the 
sheath and NG (possibly including the positive column). 

At the end of the cathode sheath it must be verified the 
following boundary condition by the displacement field D 



n.(Dp-D,) = cr. 



(29) 



Here, a is the surface charge density accumulated at the 
boundary surface and n is the normal to the surface. In 
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more explicit form, 



£pEp — EsEs — a. 



(30) 



Here, and Sp are, resp., the electrical permittivity of 
the sheath and the positive column. We have to solve 
the following algebraic system of equations 



£pEp ~ EgEg — a. 
They give the electric field strength in each region 



E., 



Ep = 



Va 

L 



la 



(31) 



(32) 



Here, we define a = 1 



Recall that in DC 



case, £p = 1 — and Sg — Sq, with uip denoting the 
plasma frequency and Ven the electron-neutral collision 
frequency. In fact, our assumption Es — £q is plainly jus- 
tified, since experiments have shown the occurrence of 
a significant gas heating and a corresponding gas den- 
sity reduction in the cathode fall region, mainly due to 
symmetric charge exchanges processes which lead to an 
efficient conversion of electrical energy to heavy-particle 
kinetic energy and thus to heating 

Two extreme cases can be considered: i) ujp > Vem 
implying Sp < 0, meaning that TcoU > Tpiasma, i-e, non- 
collisional regime prevails; ii) LUp < Ven, £p > 0, and then 
TcoU > Tpiasma, i-S, coUisional regime dominates. 

From the above Egs. [5^we estimate the field inversion 
should occurs for the condition 1 — ^ = 0, which give 
the position on the axis where field inversion occurs: 



In 



UJ. 



en 
~2' 



(33) 



From Eq. [33] we can resume a criteria for field rever- 
sal: it only occurs in the non-collisional regime; by the 
contrary, in the coUisional regime and to the extent of 
validity of this simple model, no field reversal will occur, 
since the slow electrons scattering time inside the well 
is higher than the the well lifetime, and collisions (in 
particular, coulombian collisions) and trapping become 
cornpetitive processes. A similar condition was obtained 
in |27| when studying the effect of electron trapping in 
ion-wave instability. Likewise, a self-consistent analytic 
model [5I have shown that at at sufficiently high pressure, 
field reversal is absent. 

Due to the accumulation of slow electrons after a dis- 
tance = L — Iq, real charges accumulated on a sur- 
face separating the cathode fall region from the negative 
glow. Naturally, it appears polarization charges on each 
side of this surface and a double layer is created with 



a surface charge - 
on the anode side. 
e = eo(l + Xe), X< 



-a[ < on the cathode side and 
But, cr' = (P • n), P = eoXeE with 
denoting the dimensionless quantity 



called electric susceptibility. As the electric displacement 
is the same everywhere, we have Dq = Di — D2. Thus, 
the residual (true) surface charge in between is given by 

cr = -cr;+ (72- (34) 

After a straightforward but lengthy algebraic operation 
we obtain 



B 

= SpVa^, 



where 



A = L 



£0 



and 



B 



+ 1 



£oi£s — £p) 

£s£p 



We can verify that a must be equal to 



(35) 



(36) 



(37) 



(38) 



Considering that a = e^XeE, we determine the minimum 
value of the electric field at the reversal point: 

E^ = S^. (39) 



^p Vg 



Here, Xe ~ £rw — 1, with Erw designating the relative 
permittivity of the plasma trapped in the well. From the 
above equation we can obtain a more practical expression 
for the electrical field at its minimum strength 



Em — 



T V 



(40) 



The magnitude of the minimum electric field depends 
on the length of the negative glow Iq- This also means 
that without NG there is no place for field reversal, and 
also the bigger the length the minor the electric field. 
The length of the negative glow can be estimated by the 
free path length Iq of the fastest electrons possessing an 
energy equal to the cathode potential fall value eVa- 

7o (iVF(w)) 

Here, w is the electrons kinetic energy and NF{w) is 
the stopping power. For example, for He, it is estimated 
plo = 0.02eK 01 (in cm.Torr units, with K in Volt). 
We denote by Uew the density of trapped electrons and 
by Tew their respective temperature. Altogether, Uep and 
Tep are, resp., the electron density and electron temper- 
ature in the negative glow region. 

By other side, we can estimate the true surface charge 
density accumulated on the interface of the two regions 
by the expression 



"=I = 



(42) 



Here, Q is the total charge over the cross sectional area 
where the current fiows and is the width of the po- 
tential well. 
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TABLE L Data used for E/p = 100 V/cm/Torr. 
Cross sections and electron temperatures are taken 
from Siglo Data base, CPAT and Kinema Software, 
http://www.Siglo-Kmema.com 



Gas 


Te (eV) 


a (10-'" cm^) 


Ar 


8 


4.0 


He 


35 


2.0 


O2 


6 


4.5 


N2 


4 


9.0 


H2 


8 


6.0 



A. Instability and width of the potential well 

From Eqs.[38land l42lit is easily obtained the trapping 
well width 

Ae = - , f\ ■ (43) 

It is expected that the potential trough should have a 
characteristic width of t he ord er in between the electron 
Debye length {Xoe — ^J^^^) ^'^d the mean scattering 
length. Using Eq.[43l in a He plasma and assuming Va = 
1 kV, = 1 m and Uen = 1.85 x 10^ s'^ (with Te = 0.03 



eV) at 1 Torr [n = 3.22 x 10^'^ cm"^) wc estimate « 
2.6 X 10-3 cm, while the Debye length is Xoe = 2.4 x 10"^ 
cm. So, our Eq. 33] gives a good order of magnitude for 
the potential width, which is expected to be in fact of 
the same order of magnitude than the Debye length. 

Table I present the set of parameters used to obtain 
our estimations. We give in Table II the estimate of the 
minimum electric field attained inside the well. The first 
field reversal at ~ iNn corresponds to the maximum 
density Uew ^ 'T-ep [lOjllMl- So, the assumed values for 
the ratio of electron temperatures and densities of the 
trapped electrons and electrons on the NG are typical 
estimates. 



TABLE 11: Minimum electric field at reversal point and well 



width. Conditions: He gas, 
kV, = 0.1, ^ = 10. 



1 Torr, lo = 20 cm, Va 



Em (V.cm-^) 


(cm) 


< -2.5 


2.6 X 10-^ 



It can be shown that there is no finite configuration 
of fields and plasma that can be in equilibrium without 
some external stress [2^. Consequently, this trough is 
dim to be unstable and burst electrons periodically (or 
in a chaotic process), releasing the trapped electrons to 
the main plasma. This phenomena produces local pertur- 
bation in the ionization rate and the electric field giving 
rise to ionization waves (striations). In the next section, 
we will calculate the time of trapping with a simple Brow- 
nian model. 

From Eq. [33] we calculate the cathode fall length for 
some gases. For this purpose we took He and H2 data 
as reference for atomic and molecular gases, resp. The 
orders of magnitude are the same, with the exception of 
Ar. Due to Ramsauer effect direct comparison is difficult. 

In Table HI it is shown a comparison of the experimen- 
tal cathode fall distances to the theoretical prediction, as 
given by Eq. 1431 Taking into account the limitations of 
this model these estimates are well consistent with ex- 
perimental data [30| . 



TABLE 111: Comparison between theoretical and experimen- 
tal cathode fall distance at p=l Torr, _E/p=100 V/cm/Torr. 
Experimental data are collected from Ref. 



Gas 


er (cm) 


C^' (cm) 


Ar 


7.40 


0.29 (Al) 


He 


1.32 


1.32 (Al) 


H2 


0.80 


0.80 (Cu) 


N2 


0.45 


0.31 (Al) 


Ne 


0.80 


0.64 (Al) 


O2 


0.30 


0.24 (Al) 



B. Lifetime of a slow electron in the potential well 

The trapped electrons most probably diffuse inside 
the well with a characteristic time much shorter than 
the lifetime of the through. Trapping can be avoided 
by Coulomb collisions [23| or by the ion-wave instabil- 
ity, both probably one outcome of the stress energy un- 
balance as previously mentioned. We consider a simple 
Brownian motion model for the slow electrons to obtain 
the scattering time r, and the lifetime T of the well. A 
Fermi-like model will allow us to obtain the slow electron 
energy distribution function. 

Considering the slow electron jiggling within the well, 
the estimated scattering time is 



(AO: 

-De 



(44) 
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Here, is the electron diffusion coefficient at thermal 
velocities. 

The fluctuations arising in the plasma are due to the 
breaking of the well and we can estimate the amplitude 
of the fluctuating field by means of Eq. |40l We obtain 



SE,-, 



el 



Then, we have 



Em Iq 



(45) 



(46) 



In Table IV we summarize scattering and trapping 
times for a few gases. 



C. Power-law slow electrons distribution function 

As slow electrons are trapped by the electric field inver- 
sion, some process must be at work to pull them out from 
the well. We suggest that fluctuations of the electric field 
in the plasma (with order of magnitude of £^c)act over 
electrons giving energy to the slow ones, which collide 
with those irregularities as with heavy particles. From 
this mechanism it results a gain of energy as well a loss. 
This model was first advanced by E. Fermi ^2611 when de- 
veloping a theory of the origin of cosmic radiation. We 
shall focus here on the rate at which energy is acquired. 

The average energy gain per collision by the trapped 
electrons (in order of magnitude) is given by 



Aw = Uw{t), 



(47) 



with U = £^ and where w is their kinetic energy. After 
N collisions the electrons energy will be 



'Ut 

w(t) = et exp ( — 



(48) 



with Ef being their thermal energy, typical of slow elec- 
trons. The time between scattering collisions is t. As- 
suming a Poisson distribution P[t) for electrons escaping 
from the trapping, then we state 



P{t) = exp{-t/T)dt/T. 



(49) 



The probability distribution of the energy gained is a 
function of one random variable (the energy), such as 



fwiw)dw ~ P{w < w < w + dw}. 



(50) 



This density fw{w) can be determined in terms of the 
density P(t). Denoting by ti = T the real root of the 
equation w — w{ti = T), then it can be readily shown 
that slow electrons obey in fact to the following power- 
law distribution function 



(51) 




w(eV) 



FIG. 4: Slow electrons distribution function vs. energy, for 
the same conditions as presented in Table 4. Solid curve: Ar, 
broken curve: N2. 



Like many man made and naturally occurring phenom- 
ena (e.g., earthquakes magnitude, distribution of in- 
come), it is expected the trapped electron distribution 
function to be a power-law (see Eq. [51]) . hence 1 -|- -^^f = 
n, with n = 2 ^ 4 as a reasonable guess. Hence, we 
estimate the trapping time to be 



T 

E'in 



(52) 



Fig.l shows the slow electrons distribution function 
pumped out from the well for two cases: Ar (solid curve), 
and N2 (broken curve). It was chosen a power exponent 
n = 2. Those distributions show that the higher confin- 
ing time is associated with less slow electrons present in 
the well. When the width of the well increases (from solid 
to broken curve) the scattering time become longer, and 
as well the confining time, due to a decrease of the relative 
number of slow electrons per given energy. This mecha- 
nism of pumping out of slow (trapped) electrons from the 
well can possibly explains the generation of electrostatic 
plasma instabilities. 

Note that the trapping time is, in fact, proportional to 
the length of the NG and inversely proportional to the 
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TABLE IV: Scattering time and trapping time in the well. The parameters are: E/N = 100 Td, Tg = 300 K, Va = 1 kV and 
^0 = 0.1 m. 



Gas (cm^s'^)" i^en^,(s~^) J* AC(cm) r (s) T (s) 

"Ar 2.52 X 10'' 8.10 x lO'-* 1.34 x 10"^ T?\0~kW^ 3.97 x lO"'' 

He 5.99 X 10® 2.39 x 10^ 1.54 x 10"^ 3.95 x 10"" 1.70 x 10"'^ 

N2 6.11 X 10^ 6.15 X lO'' 2.32 x 10"^ 8.81 x 10"^^ 1.64 x 10"" 

CO2 1.70 X 10" 3.60 X lO'' 6.78 x 10"^ 2.70 x 10~" 5.90 x 10"^ 



"Data obtained through resolution of the homogeneous electron 
Boltzmann equation with two term expansion of the distribution 
function in spherical harmonics, M. J. Pinheiro and J. Loureiro, J. 
Phys. D.: Appl. Phys. 35 1 (2002) 

''Same remark as in " 



electrons diffusion coefficient at thermal energies: 



VIII. CONCLUSION 



T 



(53) 



Tlie survival frequency of trapped electrons is — l/T. 
As the electrons diffusion coefficient are typically higher 
in atomic gases, it is natural to expect plasma instabil- 
ities and waves with higher frequencies in atomic gases. 
This result is in agreement with a kinetic analysis of in- 
stabilities in microwave discharges [3l|. In addition, the 
length of the NG will influence the magnitude of the fre- 
quencies registered by the instabilities, since wavelengths 
have more or less space to build-up. Table HVl summarizes 
the previous results for some atomic and molecular gases. 
The transport parameters used therefor where calculated 
by solving the electron Boltzmann equation, under the 
two-term approximation, in a steady-state Townsend dis- 
charge [3^ . 



This paper reviews, and extends when necessary, the 
physics of the electric field reversals in a glow discharge. 
We listed and studied how the structure of a glow dis- 
charge is related to the electrons and ions kinetics, de- 
scribed related experimental results, and have shown 
how the field reversals help to maintain the discharge. 
To this end different methods to describe the different 
groups of electrons generated in a glow discharge have 
been presented: the nonlocal approach of Kolobov and 
Tsendin the simple analytic fluid model by Boeuf 
and Pitchford [l3|; the theoretical analysis developed 
by Kudryavtsev and Toinova [20] , and the dielectric-like 
model [221. 
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